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SECTION  I 


INTRODUCTION 


'  In  two  previous  reports  (Schneck,  1976,  1979)  a  mathematical  theory  was 
developed  in  order  to  calculate  the  aerodynamic  loading  to  which  a  pilot  Is 
exposed  during  high  speed  ejections.  Neglecting  the  effects  of  flow  sepa¬ 
ration,  preliminary  results  suggested  that  the  pilot's  musculo -skeletal 
resistance  is  not  likely  to  withstand  the  tendency  for  dislodgement  from  a 
restraining  surface  if  he  is  ejecting  at  Mach  numbers  exceeding  0 . 72 .—.More 
specifically,  the  aerodynamic  dislodging  (windblast)  force,  F  (in  lbs.)],  that 
is  generated  at  an  altitude  of  10,000  feet  can  be  described  by  the  equation: 


1148m2 sin2a. 


(1) 


where  a  is^the  angle  at  which  a  pilot's  limb  intercepts  an  air  stream  moving 
at  Mach  ^M". 


^“Recently  (Schneck,  1979)  the  mathematical  theory  of  windblast  was  modi¬ 
fied  to  include  some  effects  of  flow  separation.  In  the  report  that  follows, 
these  effects  are  examined  in  greater  detail,  particularly  as  they  affect  the 
time-course  of  limb  dislodging  forces  after  the  onset  of  windblast. 


SECTION  II 

VORTEX  PATTERNS  FOLLOWING  DISLODGEMENT 


Consider  the  situation  depicted  in  Figure  1.  This  illustrates  the 
instantaneous  streamline  pattern  that  characterizes  the  cross-flow  over  two 
circular-cylindrical  body  segments  in  line  contact  with  one  another  at  time  t 

*  0.  That  is,  at  the  moment  this  limb-to-limb  or  limb-to-res training  surface 

configuration  intercepts  a  high-speed  cross-flow,  the  streamline  pattern  is 
initially  inviscid,  a  boundary  layer  has  had  no  time  to  gr<-  4  the  flow 
starts  out  essentially  attached.  However,  this  situation  w  .y  last 

momentarily,  since  stagnation  of  the  flow  along  the  line  contact  2 tween  the 
two  cylinders  will  result  in  the  generation  of  relatively  large  forces  which 
tend  to  pull  the  cylinders  apart.  In  fact,  equation  (1)  predicts  that,  for  a 

*  90°,  beyond  Mach  0.72  the  forces  so  generated  will  exceed  the  amount  by 
which  an  average  pilot  can  resist  the  tendency  to  have  his  limbs  pulled 

apart.  Thus,  it  is  likely  that  for  t  *  0+,  a  fraction  of  a  second  after  the 
flow  is  intercepted,  the  situation  depicted  in  Figure  2A  prevails.  Here,  the 
cylinders  are  shown  separated  by  a  very  small  distance,  b,  simulating  the 
beginning  of  flail  resulting  from  windblast. 

The  configuration  illustrated  in  Figure  2A  produces  a  venturi  effect. 

The  fluid  passing  through  the  throat  between  the  cylinders  separated  a  slight 
distance  b  is  accelerated  and  exits  the  gap  with  some  velocity,  u^.  On  the 


other  hand,  the  fluid  diverted  to  the  outside,  not  passing  through  the  throat, 
circles  around  the  cylinders  and  diffuses  sharply  on  the  downstream  side, 
decelerating  to  some  velocity  u^  «  u^.  Thus,  in  the  narrow  region,  y,  where 

a  streamline  having  velocity  u^  meets  a  streamline  having  velocity  u2,  there 

is  a  large  velocity  gradient  which  causes  the  fluid  to  roll  up  into  a  starting 
vortex,  as  shown  in  Figure  2B. 

The  starting  vortex  thus  generated  on  each  cylinder  grows  until  it  is 
eventually  shed,  leaving  behind  an  equal  and  opposite  circulation,  T,  about 
the  cylinder  as  shown  in  Figure  2C.  Note,  however,  that  this  circulation 
acts  to  retard  fluid  flow  through  the  throat  region  while  increasing  the 
flow  around  the  outer  sides  of  the  cylinders.  It  thus  works  against  the 
venturi  effect  depicted  in  Figure  2A.  The  result  is  that,  in  time,  the 
effect  of  the  circulation  is  to  make  the  fluid  velocity  on  the  outer  sides  of 
the  cylinders  greater  than  that  between  the  cylinders.  This  reverses  the 
velocity  gradient  illustrated  in  Figure  2A,  such  that  vortices  of  opposite 
spin  now  roll  up  on  the  outer  rim  of  the  cylinders  as  the  faster  moving  fluid 
encounters  the  slower  moving  streamlines  diffusing  out  on  the  downstream  side 
of  the  throat.  This  situation  is  illustrated  in  Figure  2D. 

In  turn,  the  second  set  of  vortices  generated  as  described  above  grows 
until  it  is  shed,  leaving  behind  an  equal  and  opposite  circulation  which 
effectively  cancels  the  circulation  produced  by  the  initial  set  of  vortices. 
Thus,  the  picture  essentially  reverts  back  to  that  shown  in  Figure  2A  and  the 
whole  process  starts  all  over  again.  Over  a  period  of  time,  then,  when  a 
steady-state  has  been  reached,  the  wake  behind  the  slightly  separated  double 
cylinder  configuration  can  be  described  as  containing  a  vortex  street  con¬ 
sisting  of  pairs  of  vortices  shed  periodically,  180°  out  of  phase,  and 
alternating  in  direction  of  rotation  as  illustrated  in  Figure  2E.  The 
spacing  between  rows  of  vortices  is  designated  as  shown  (h  between  rows  1  and 
2  and  3  and  4,  respectively  and  k  between  rows  2  and  3),  while  the  distance 
between  vortices  in  a  row  is  given  by  i. 

The  significance  of  the  entire  preceding  discussion  rests  in  the  reali¬ 
zation  that  the  oscillatory  nature  of  the  flow  described  results  in  a  corres¬ 
ponding  oscillation  of  the  forces  tending  to  cause  lateral  motion  of  the  two 
cylinders.  That  is,  as  vortices  in  rows  2  and  3  are  shed,  the  circulation 
left  behind  causes  u2  to  fall  well  below  u^  so  the  corresponding  pressure  p2 

>>  p.^  and  there  is  a  tendency  for  the  cylinders  to  be  driven  apart  from  one 

another.  Conversely,  as  vortices  in  rows  1  and  4  are  shed,  the  resulting 
circulation  leads  to  u^  >>  u^,  hence  p2  <<  and  the  cylinders  are  driven 

towards  one  another.  The  net  effect  of  the  double  vortex  street  in  the  wake 
is  to  produce  a  lateral  vibration  of  the  cylinders  once  they  have  separated 
due  to  flow  stagnation  forces.  In  the  section  that  follows  the  theory  for 
calculating  the  oscillating  forces  is  developed  following  the  method  outlined 
by  Landweber  (1942)  and  Milne-Thomson  (1955). 
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Figure  2:  The  Establishment  of  a  Double  Vortex  Street  Pattern  Following 
the  Interaction  of  the  Configuration  Depicted  in  Figure  1  with 
a  High-Speed  Cross-Flow. 


SECTION  III 


DEVELOPMENT  OF  THEORY  FOR  CALCULATING 
OSCILLATING  FORCES  PRODUCING  LATERAL  DISPLACEMENT 


Referring  to  Figure  3,  consider  a  circular  cylinder  of  radius  "a"  lying 
perpendicular  to  a  uniform  stream  of  velocity  V  and  exposed  to  a  circulation, 
T,  of  strength  k,  such  that  T  •  2nx.  The  lef t-to-right  streaming  motion  past 
the  cylinder  is  given  by  the  complex  potential: 


v1  -  -V[z  +  —3  ,  z  -  x  +  iy  -  re10  (2) 

while  the  clockwise  circulation  of  strength  <  about  the  cylinder  is  given  by: 

w2  -  -itc  log  j.  (3) 

Combining  these  motions  yields  the  complex  potential: 


2 

w  ■  w  +  w„  *  -V[z  +  — ]  -  Ik  log  — .  (4) 

1  2  z  a 

Let  there  be  a  force  f  ■  F  +  iF  acting  per  unit  length  on  the  cylinder 

x  y 

as  a  result  of  the  pressure  distribution  generated  by  the  flow  defined  in 
equation  (4) .  Then: 


F 


•fa  -f, 


p[idx  +  dy] 


(p)d[x  -  iy]  *  fipdz 


Now,  the  steady-state  form  of  Bernolli's  equation  is, 

1  2 


p  +  j  pq 


constant,  C 


(5) 

(6) 


where : 


2 


q 


dw  dw 


But,  since  w 


+  ii^  and  w  *  $  -  iiji,  and  since  the 


surface  of  the  cylinder  is  the  streamline  ■  0,  we  note  that,  on  the  surface 
w  *  w,  and  so  dw  *  dw.  Moreover,  if  interest  is  focused  only  in  the  pres¬ 
sure,  p,  that  is  due  to  the  flow,  q,  i.e.,  we  examine  from  the  reference 
point  p  (or  p  minus  some  stagnation  pressure  p  )  ■  0  when  q  *  0,  then  C  ■  0, 
and  equation  (5)  becomes: 


; ./«-  i  *>« ./«-  ;>  if  « 1 


From  equation  (4) : 


(7) 
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Figure  3:  Circular  Cylinder  of  Radius  "a".  Lying  Perpendicular  to  A 
Uniform  Stream  of  Velocity  V,  and  Exposed  to  a  Circulation 
T,  of  Strength  k,  such  that  T  =  2mc. 


U  ■  - {v[i  -  +  u 

z 


2 

Thus,  neglecting  compressibility  effects,  ?  ■  -  {V  -  y  +  ^-}2dz  (8) 

z* 

Equation  (8)  is  a  complex  integral  to  be  taken  around  the  contour  of  the 
cylinder.  The  integrand  can  be  expanded  to  give: 


v2  +  liV< 
z 


2aV  k2  2ia2V<  4”2 


2  2 
z  z 


3  4 

z  z 


which  is  a  specific  form  of  the  more  generalized  power  series  expansion: 
00 

f(z)  -  ^  aQ(z  -  zo)n,  where. 


n--<» 


(?-z  > 
o 


“T-  In  this  case,  z  *  0,  and  -4  <  n  <  0,  with: 
n+1  o  —  — 


a  *  V 
o 


*-l  *  2^i 


f(C)dC  ■  sum  of  residues  inside  the  contour  ■  2iV< 


2  2  2 
a_2  *  -(2aV  +  O 

a_2  *  -2ia2Vtc  and, 
4„2 

a_4  -  a  v 


'•i 


Thus,  pf(z)dz  *  2iria_^  ■  27ri(2iVic),  and  so: 


f  -  -  ±1  (2iri)(2iVic)  -  2iprVie 

or,  ?  »  F  +  iF  *0  +  i(2pirVK).  Finally,  then: 
x  y 

F  *  0,  and 
x 

f  -  rpv 
y 


(9) 
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Note  the  absence  of  a  drag  force,  F  ,  because  the  flow  Is  Initially  con¬ 
sidered  to  be  irrotational .  Furthermore,  note  that  the  "lift"  force,  F^, 

depends  on  the  circulation,  T,  about  the  cylinder.  This  circulation  may  be 
described  in  terms  of  the  vortex-shedding  characteristics  of  the  double¬ 
cylinder  configuration  as  follows:  consider  a  point  P  located  at  coordinate 

z  ■  x  +  iy  as  shown  in  Figure  4  relative  to  the  Row  1  vortices  shed  as 
o  o  Jo 

illustrated  in  Figure  2.  The  velocity  at  point  P  due  to  any  vortex,  say  Q, 
in  row  1,  may  be  obtained  from  the  complex  potential  function  for  that  vortex 
filament  of  strength 


w  *  ”i,ci  l°g(z  “  ZQ)  (10) 

4Q 

where  (z  -  z  )  *  re  defines  the  location  (z)  of  the  vortex  relative  to 
o 

arbitrary  point  P(xQ,yQ)  in  the  x-y  plane.  The  negative  sign  in  equation 

(10)  indicates  that  the  circulation  around  the  filament  is  clockwise  as  shown 
in  the  Figures.  From  equation  (10),  we  may  write: 


w  ■  -i<^  log(rei9)  *  -iK^log  r  +  i6]  -  +  lip. 

Thus,  <|>  -  k^0  and  ip  »  -k^  log  r. 


from  which:  v  =-|£»0=*-||--^ 
r  dr  oo  r 

i  d&  _  =  ii 

V8  3  '  r  30  “  r  *  3r 


(11) 


and  so  it  is  observed  that  the  direction  of  motion  at  P  due  to  the  vortex  at 
Q  is  at  right  angles  to  the  line  QP  or  r  (since  vr  *  0) ,  with  speed  K^/r,  in 

the  sense  given  by  the  rotation  of  the  vortex  at  Q  (hence  the  negative  sign 
in  the  equation  for  vQ) .  We  may  rewrite  the  equation  for  vQ  as: 

v@  "  VP  "  "  r  V2*  (12) 


where  »  2jtk^  is  the  circulation  of  each  vortex  in  row  1.  If  we  now 

examine  a  "mirror"  vortex  located  at  the  corresponding  point  Q'  in  row  1,  it 
is  noted  that  its  effect  at  P  is  to  induce  a  y-velocity  component  that  can¬ 
cels  the  one  induced  by  the  vortex  at  Q,  and  an  x-velocity  component  that 
doubles  the  one  already  there  due  to  the  vortex  at  Q.  Hence,  the  net  velocity 
at  P  due  to  vortices  at  Q  and  Q'  is  in  the  negative  x-direction  and  has 
magnitude: 
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Velocity  Induced  at  an  Arbitrary  Point,  P,  due  to  the  Vortex 
Stream  Generated  in  Row  1  Isolated  from  Figure  2. 


X 


'57  7  (2>  8109 


1  -  rsin6 

"  2ir  Z  2 
r 


_  _ 


y  + 


9£ 


2  * 


Generalizing  the  above  discussion  to  any  n'th  pair  o£  vortices  in  row  1,  we 
2  l  2  2 

may  write:  r  ■  (n£  -  -j )  +  y  ,  and  summing  the  effects  of  all  the  vortices 

in  the  row,  we  conclude  that  the  velocity  induced  at  P  by  the  vortex  stream 
in  row  1  is: 

r,  j.  r. 


£1=1 


(ni,  -  -|)2  +  y2 


-E4 

n*l  it l 


(» 


r 

_i  iz. 

2£  Ttl 


(2n  -  l)2  +  ^ 
£ 


r. 

— i.  tanh  £ 
21  cann  1 


or, 


(13) 


where 


*  ir£ 


(2n  -  l)2+^_ 


tanh(y 


(Gradshteyn  and  Ryzhik,  1965,  page  36).  Performing  a  similar  analysis  for  a 
point  P'  located  as  shown  in  Figure  4,  one  obtains  (Landweber,  1942): 


V  *  "  2£  COth  l 


(14) 


Using  equations  (13)  and  (14),  one  can  now  write  relationships  for  the 
velocity  of  vortices  in  any  row  of  the  street  shown  in  Figure  2E.  That  is, 
vortex  in  row  4  is  travelling  with  a  speed  induced  by  the  vortices  in  rows  1 
2,  and  3.  Hence: 


^2  .  ,  irh  .  ^2  ,  (h-Hc)iv 


-  ii  coth  sum. 
2  £  £ 


(15) 


where  ^  t^e  circulation  of  each  vortex  in  rows  2  and  3.  Also,  a  vortex 
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in  row  3,  is  travelling  with  a  speed  induced  by  the  vortices  in  rows  1,  2, 
and  4.  Hence: 


r  r 

1 1  u  irh  ,  2  . .  irk 

'  n  tanh  T  +  IT  coth  T 


F1  .  .  ir(k+h) 

'  21  tanh  “l 


(16) 


Examination  of  equations  (15)  and  (16)  clearly  reveals  that  F^  and 

could  be  calculated  if  we  could  say  something  about  v^,  v^,  v^,  v^,  l,  h, 

and  k.  From  stability  considerations  we  do  know  that  *  v2  ”  v3  “  v4»  but 

more  information  must  be  obtained.  This  can  be  achieved  most  conveniently  by 
performing  a  series  of  relevant  and  carefully  designed  experiments.  That  is 
the  subject  of  the  section  which  follows.  But  before  proceeding,  it  is  of 
Interest  to  use  some  related  experimental  data  acquired  by  Landweber  quite  a 
long  time  ago.  From  his  studies,  Landweber  (1942)  suggested  that  the  following 
are  reasonable  values  for  the  parameters  in  equations  (15)  and  (16) :  for 
small  b,  (  0  <  b  <  .125a),  and  for  steady  turbulent  flow, 

|  «  1  ,  l  -  rrh  ,  v2  *  v2  *  V3  ’  V4  =  4 

and  h  *  (1.3) 2a. 

Substituting  into  (15) : 

v  '  ,(1.3)  (2)  (2a)  [-r2  +  f2  ta"h(2>  '  fl  C°th(3)1 

thus,  5 . 2irav  =*  -.7616r2  +  .9640r2  -  1.00501^. 

Substituting  into  (16): 

v  -  27r(-£13-y2-  [-rx  tanh(l)  +  r2  coth(l)  -  rx  tanh(2)] 
thus,  5 . 2irav  -  -.7616T  +  1.3130r2  -  .9640F;L. 

Rewriting: 

.2024r2  -  l.OOSOrj^  -  1.3iraV 
and,  1.3130r2  -  1.7256ri  -  1.3iraV 
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Solving  simultaneously:  .lOlAr^  -  ■  1.2935iraV 

.7609r.  -  r,  -  .7534iraV 

_ _ 

Subtract:  -.5595^  a  .5401iraV 

Thus,  r,  a  -3. 0327aV, 

(17) 

and,  rL  ■  -4.6745aV, 

where,  in  the  steady  state,  ^  is  the  maximum  circulation  in  rows  2  and  3 

driving  the  cylinders  apart  and  F^  is  the  maximum  circulation  in  rows  1  and  4 

driving  the  cylinders  together.  The  negative  sign  is  appropriate  for  rows  1 
and  3,  while  it  changes  to  a  positive  sign  for  rows  2  and  4. 

SECTION  IV 

CALCULATION  OF  OSCILLATING  LATERAL  FORCES 


At  this  point  we  have  some  estimate  of  the  information  necessary  to 
evaluate  equation  (9)  in  order  to  determine  the  forces  which  cause  the 
cylinders  to  oscillate  after  they  have  separated  a  small  amount.  To  be 
consistent  with  an  earlier  report  (Schneck,  1979)  it  is  assumed  that  the 
radius  of  the  cylinder  is  1.5  inches  (characteristic  of  the  human  forearm) 
and  it  has  a  length  of  one  foot.  Furthermore,  an  ejection  taking  place  at  an 
altitude  of  10,000  feet  is  considered,  where 

*  0.7385,  p  a  0.002378  slug/ft^  and  c,  the  local  speed  of  sound  is 

Po  ° 

1,078  ft/sec.  Thus, 

F  *  r(£-)p  -  c  *  1.8931  r  M  sin  a  (18) 

v  'p  "o  c 
o 

where:  V  *  Uq  sin  a  =*  magnitude  of  the  cross-flow  velocity, 

Uq  =*  magnitude  of  free-stream  air  flow  velocity, 

a  =  angle  of  attack  between  cylinder  axis  and  direction  of 
air  flow  velocity,  and, 

M  *  U  /c  ■  Mach  number  of  the  flow, 
o 

Substituting  f ^  from  equation  (17)  into  equation  (18),  with  a  ■  1.5 


inches,  we  obtain  the  maximum  force  driving  the  cylinders  apart,  which  is: 


1 


.  ,  U  sina 

F2  -  1.893ir2Msina  =■  1.8931(3. 0327)^  S~- —  (1078)Msina,  or, 

F2  *  773.6  M2sin2a  (19) 

Similarly,  substituting  F^  from  equation  (17)  into  equation  (18)  we  obtain 

the  maximum  force  driving  the  cylinders  together,  which  is: 

-  .  U  sina 

Fx  -  1.8931[4.67A5(-^)  9~g—  (1078)Msina] ,  or, 

Fx  -  1192.4  M2sin2a  (20) 


The  time-course  of  aerodynamic  forces  that  prevail  some  time  following 
the  onset  of  windblast  thus  oscillates  between  the  extremes  defined  by  equa¬ 
tions  (19)  and  (20) .  The  frequency,  m,  of  such  oscillation  depends  on  the 
vortex  shedding  frequency,  which  is  usually  expressed  in  terms  of  the  dimen¬ 
sionless  Strouhal  number: 


c  ,  °(2a) 
S  V 


(21) 


Thus,  m  *  -jg.  But,  this  vortex  shedding  frequency  can  also  be  written  in 


terms  of  the  number  of  times  per  second  that  vortices  in  a  given  row,  j ,  pass 
any  particular  point,  i.e.. 


V  -  v 

»  "  ••£  2  •  (22) 

V  Uo 

Writing,  as  before,  v.  ;  -r,  V  *  U  sin  a,  l  *  irh,  h  «  2.6a,  M  *  — ,  a  -  1.5 

j  4  o  c 

inches,  and  c  ”  1078  ft/sec,  equation  (22)  reduces  to: 

V  -  g  UQsin  a  M(1078)sin  a 
”  ’  —  ‘  n(2.6)a  '  ,(2.6)%5i.  •  ot- 


m  ■  791.86  M  sin  a  per  second  (23) 

with  a  corresponding  Strouhal  number: 

|  V(2a) 

S  ■  ~(y -6~aV  -  .1836  characterizing  the  flow. 


Now  write:  u  *  2irm  ■  circular  frequency 


F  +  F 
2  1 

A  ■  - 2 - ■  amplitude  of  the  oscillation 

A  ■  F„  -  A  ■  A  -  F,  ■  mean  of  the  oscillation, 
o  2  l 

Then,  from  equations  (19),  (20)  and  (23): 

id  ■  4975.38  M  sin  a, 

A  -  983  M2 sin2a,  and 

A  -  -209.4  M2 sin2a 
o 

so  that  the  oscillating  forces  generated  by  the  double  vortex  street  described 
earlier  can  be  defined  by  the  equations: 

t  >  0+  F(t)  "A  +  Acoswt 

( s teddy  ^ 

-  M2sin2a[-209.4  +  983cos(4975.4Msina)t]  (24) 

t  =  0  F(0)  -  1148M2sin2a  (25) 

where  the  negative  sign  corresponds  to  the  cylinders  being  driven  together 
and  the  positive  sign  represents  the  cylinders  being  driven  apart. 

Now,  it  is  reasonable  to  expect  that  the  transition  from  equation  (25) 

at  t  •  0  to  equation  (24)  in  the  steady  state  for  t  >  0+  does  not  occur 
instantaneously  or  discontinuously,  but  rather,  is  characterized  by  some  time 
constant,  t.  Due  to  the  physics  of  the  process,  however,  one  can  surmise 
further  that  x  is  very  small,  perhaps  on  the  order  of  one-quarter  or  less  of 
the  period  defined  by  the  reciprocal  of  equation  (23) .  Such  a  small  time 
constant  would  cause  the  transient  stage  of  the  flow  adjustment  to  last 
effectively  one  complete  cycle,  since  after  a  lapse  of  time  equal  to  four 
time  constants  the  transient  components  of  the  process  will  have  decayed  to 

1/4 

some  27.  of  their  initial  value.  Thus,  let  x  *  -;v, — -z —  and  combine 

791.86Msina 

equations  (24)  and  (25)  in  the  form: 

_  t 

F(t)  *  M2sin2a{374.4e  T  -  209.4  +  983  cos [ 4975. 4Msina]t}  (26) 

observe  that  equation  (26)  reduces  to  equation  (25)  for  t  *  0,  and  equation 
(26)  reduces  to  equation  (24)  as  t  +  ®.  In  a  practical  sense,  however,  the 
leading  term  of  equation  (26)  is  down  to  2%  of  its  initial  value  by  the  time 
t  *  4x,  which  is  after  one  complete  cycle  has  transpired. 

Figure  5  shows  equation  (26)  plotted  with  M  sina,  acting  as  a  parameter. 
Essentially,  M  sina  may  be  redefined  by  some  new  variable,  say  M^,  repre- 
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senting  the  cross-flow  Mach  number  of  the  air  passing  over  the  limbs  of  the 
ejection  seat  occupant.  Note  that  Mg  may  vary  over  different  portions  of  the 

body  that  intercept  the  flow  at  different  angles  a. 

Equation  (25),  valid  at  t  ■  0,  (or,  the  limit  of  equation  (26)  as  t 
0)  plots  along  the  vertical  axis  at  the  origin  of  the  coordinate  system,  and 
defines  the  limb-dislodging  fore*  that  prevails  for  the  situation  depicted  in 
Figure  1.  Equation  (24),  (or,  the  limit  of  equation  (26)  as  t  ®)  then 
describes  the  oscillating  lateral  forces  that  are  associated  with  the  double 
vortex  street  pattern  which  develops  after  some  time  in  accordance  with  the 
events  illustrated  in  Figure  2  for  t  >  0.  The  parameter  M  sina  is  examined 
only  for  cross-flow  Mach  numbers  M.  1,  i.e. ,  in  the  sub-and-early-transonlc 
range . 

Observe,  first  of  all,  that  the  frequency  of  the  oscillating  lateral 
forces  Increases  linearly  with  Mg  according  to  equation  (23)  (plotted  in 

Figure  6),  and  that  it  becomes  rather  substantial  at  higher  cross-flow  Mach 
numbers.  This  may  explain,  in  part,  the  observed  and  often  harmful  vibration 
which  takes  place  between  the  head/neck  of  the  pilot  and  the  ejection  seat 
head  rest  during  higher-speed  ejections.  Even  allowing  for  the  inertial 
resistance  of  the  mass  of  the  head,  it  is  likely  that  vibrations  on  the  order 

2 

of  10  cycles  per  second  may  accompany  ejections  occurring  above  Mg  -  0.7, 
where  the  forces  generated  exceed  the  limits  of  musculo-skeletal  resistance. 
Writings 


F(t)  -  M  a  (t),  (27) 

S  9 

assuming  a  simple,  unrestrained,  two  dimensional,  rectilinear  situation,  one 
may  estimate  the  kinematics  of  the  ensuing  motion  of  a  body  segment  of  mass 
Mg,  subjected  to  the  force  distribution  defined  by  equation  (26).  Thus, 

a  (t)  -  SS£l  .  {374>4e(-3167.44tMsincO  _  2Q9>4 

S  M  M 

s  s 

+  983cos[4975.  4Msinct]t)  (28) 

Integrate  equation  (28)  to  get  velocity: 


vg(t) 


u2  ,  2 
M  sin  a 

M 


1182e("3167,44tMsina) 

Msina 


209. 4t 


+  .1976 
Msina 


sin[4975.4Msina]t 


+  ci> 


(29) 
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(onoobs  aad  saioAO)  w 
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1182 

where:  *  MaTnn  ^  assumed  that  vg(0)  *  0.  Going  one  step  further, 

integrate  equation  (29)  to  get  corresponding  limb  displacement: 


Vt) 


„2  .  2  ,_««  1r.-4  (-3167.44tMsina)  , 
M  sin  a  ,.3732  x  10  e  ,.2 
- - -  { - - - - - 104. 7t 


sin^o 


.3971  x  10  . .  .  .1182t  .  „  , 

’  "ATA —  C08f4975'4M8±na]t  +  iteiSS"  +  C2} 

.3732  x  10'4  .  .3971  x  IQ"4 . .  .  . 

where:  * - 2 - z - + - ~2 - 2 -  ^  iC  is  ^urt^er  assumed  that 

“  M  sin  a  n  sin  a 


6  (0)  -  0. 

s 


Writing  t  *  —  T,  where  t 


[3.1371  x  10 
Msina 


equation  (30)  reduces  to: 


6  (t)  -  {.3732  x  10-4e  T-  .1044  x  10-4£)2 


-.3971  x  10‘4cos(1.5708>E- +  .3732  x  10~4  - 

T  T 


+  .0239  x  10~4}, 


5  (n)  “  {.3732e"n-.3971  cos^  n) 

S  M  i 

s 


-.1044n  +  .3732n  +  .0239) 

where:  n  *  Equation  (31)  defines  an  exponentially  decaying  cosine  wave 
super-imposed  on  a  parabola  represented  by  the  last  three  terms.  One  com- 

4 

plete  cycle,  (4  time  constants)  is  shown  in  Figure  7,  where  <5g(n)Mg10  is 

plotted  against  n  for  0  £  n  <_  4.  Observe  that  during  this  early  phase  of  the 
ejection  process,  the  tendency  for  the  displacement  of  a  limb-segment  is 
directed  away  from  a  restraining  surface,  i.e.,  6g(n)  is  positive  at  least  up 

to  ri  ■  3.  While  the  actual  displacements  predicted  by  equation  (31)  are 
small  (typically  on  the  order  of  0.20  mm),  our  concern  in  this  simple 
example  is  less  with  specific  numbers  than  with  physical  trends.  That  is, 
the  fact  that  <S  (n)  is  initially  positive  is  significant  because  when  the 

peripheral  limbs  (arms  and  lower  legs)  become  separated  from  the  torso  area 


in  a  situation  where  the  seat  provides  little  or  no  restraint,  drag  forces 
are  more  likely  to  cause  an  uncontrolled  flail,  rather  than  a  simple  lateral 
vibration  (such  as  might  occur  between  the  head  and  head  rest).  Moreover, 
this  theory  suggests  what  the  evidence  has  revealed  to  be  the  case,  i.e., 
that  contrary  to  what  one  might  suspect  a^  priori,  there  is  a  violent  forward 
head  rotation  (away  from  the  restraining  head  rest)  at  the  initial  instant  of 
ejection  -  with  a  strong  likelihood  of  incipient  neck  injury.  Of  course, 
musculo-skeletal  resistance  has  not  been  taken  into  account  to  this  point  in 
the  analysis;  but  Figure  5  further  reveals  the  likelihood  that  the  oscil¬ 
lating  forces  generated  beyond  *  0.70  exceed  at  least  the  average  grip 

retention  force  reported  by  Horner  and  Hawker  (1973)  for  significant  portions 
of  the  pulse  cycle.  It  thus  follows  that  flail  seems  somewhat  inevitable  at 
these  higher  Mach  numbers  and  angles  of  attack,  regardless  of  the  pilot's 
efforts  to  "hold  on."  In  any  case*  equations  (28)-(31)  are  only  simplifi¬ 
cations  intended  to  reveal  a  trend.  More  exact  kinematic  Information  can  be 
obtained  by  exercising  the  Aerospace  Medical  Research  Laboratory's  Articu¬ 
lated  Total  Body  Model  (ATB)  with  equation  (26)  providing  input  data  relative 
to  the  aerodynamic  loading  that  prevails  during  a  given  ejection  situation. 

In  fact,  even  equation  (26)  is  an  approximation  (although  realistic)  based 
on  the  results  of  Landweber  (1942).  It  is  thus  desirable  to  back  up  still 
further  to  equations  (9) ,  (15)  and  (16)  and  to  accumulate  somewhat  more 
relevant  experimental  Information  concerning  the  wake  that  develops  down¬ 
stream  of  a  double-cylinder  configuration.  To  do  so,  we  shall  exploit  the 
principle  of  the  Hydraulic  Analogy  as  described  below. 


SECTION  V 

PRINCIPLE  OF  THE  HYDRAULIC  ANALOGY 


The  idea  behind  the  hydraulic  analogy  concept  is  to  be  able  to  deduce 
what  goes  on  in  a  two-dimensional  high-speed  compressible  flow  (e.g.,  air¬ 
craft  ejections)  from  measurements  made  in  an  incompressible,  open-channel 
flow.  The  analogy  itself  is  based  on  the  correspondence  between  pressure 
waves  in  a  compressible  gas  and  gravity  waves  on  the  free  surface  of  a 
liquid  (see,  for  example,  Shapiro,  1949).  For  example,  consider  the  flow  of 
a  liquid  in  an  open  channel  of  arbitrary  cross-section  as  shown  in  Figure 
8.  Here,  the  equation  of  the  cross-section  is  taken  to  be  a  power  law  of 
the  form: 

x(y)  -  By3,  (32) 

which  is  assumed  to  be  uniform  in  the  x-direction  (coming  out  of  the  paper) . 
Let  the  maximum  water  level  at  some  arbitrary  x-section  be  corresponding 

to  a  surface  channel  half-width  g^.  Then  the  fluid  filled  cross-sectional 

area,  A^  at  this  x-station  is: 


tion 


’8*» 

ax  -  2  £cL  -  x(y)]<*y 
o  A 

•8,  f®l  ft 

-  2  1  C-jdy  -  2  1  Bypdy 

*  O  *'  O 


2S  l8i.iiz!!i|gi 

^ly,o  8+1  o 


51  8181 


2^i8i  "  e  ^e+i 
81 

2Cl8l6  +  2ClSl  -  2C1&1  2^85^ 

6+1  cb+db3 


therefore  ^  ^  ^  ^  •  (33) 

Now,  due  to  the  existence  of  8ravity  wave  motion,  the  fluid-filled  cross- 
sectional  area  at -some  other  x-station  will  be  different  from  A^,  say, 

8+1 

A  *  - 5L  c  e  (34) 

A2  8+1  b1/8  *2  *  ^  ' 


But,  for  an  incompressible  fluid,  continuity  requires  that: 

A1V1  *  A2V2  (35) 

where  and  are  the  mean  flow  velocities  at  stations  1  and  2 ,  respec¬ 
tively.  Thus, 


. 

1/8  6 


II  ,  ll  ,  _JL_  _2_  8  iiic!, 

VL  A2  8+1  b1/8  *1  8  2^2 


v,  5, 


2  ^  6 


-  *  (-i) 
V  V  ' 
1  **2 


In  order  to  have  equation  (36)  be  analo8oas  to  a  correspondin8  compressible 
flow  situation,  we  note  that  Che  continuity  equation  for  the  two-dimensional 
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steady  flow  of  a  compressible  fluid  through  a  region  of  constant  cross-sec¬ 
tion  is: 


P 31  ^2^2’  * 


(37) 


(38) 


Thus,  in  the  hydraulic  analogy,  equating  equations  (38)  and  (36)  yields: 


(^)  6  =  ^  (39) 

*2  P2 

so  that  depth  changes  due  to  wave  motion  in  an  incompressible  open  channel 
flow  are  proportional  in  the  analogy  to  density  changes  due  to  wave  motion 
in  a  compressible  gas-dynamic  flow.  Carrying  the  analogy  one  step  further, 
consider  Bernoulli’s  equation  applied  on  the  free  surface  of  the  open  channel 
flow: 


(40) 


An  "average"  form  of  equation  (40)  may  be  written  more  generally  for 
a  channel  of  arbitrary  cross-section  by  introducing  the  concept  of  the  hy¬ 
draulic  mean  depth,  5*.  This  is  essentially  a  modified  parameter  that 
accounts  for  the  fact  that  the  cross-section  of  the  open  channel  may  not  be 
rectangular.  The  hydraulic  Mean  Depth  is  defined  as  the  ratio  of  the  fluid 
filled  cross-sectional  area,  A,  at  any  x-station,  to  the  wetted  perimeter, 

P  at  that  same  station.  The  latter  is  defined  as: 
w 

pw  -  zJ*n/i  +  (^)2  dy  (41) 


In  order  to  evaluate  equation  (41),  we  first  make  some  observations  con¬ 


cerning  the  derivative  (4^0.  That  is,  let  a 


dy 


and  t 


8, 


Then, 


m 


dy  de  dy  g^ 


de  ’  so> 


5 

—  *  4—  —  ■  —  —  .  so.  7^*  ■  ~  4^-  ■  order  (— )  .  For  a  shallow  (small  C  ) 

dy  8o  d£  5  «m 

wide  (large  g  )  channel,  then,  where  —  <<  1,  one  may  certainly  conclude 
®  £*  8m 

that  (4^) 2  *  0(— )^  is  negligible  compared  with  unity.  The  neglect  of  this 
derivative  in  determining  Pw  is  thus  a  reasonable  approximation,  so  we  may 
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write,  for  a  wide,  shallow  channel 


P 

w 


From  equation  (33) ,  A  = 


finally  be  expressed  as: 


dy  ■  2y|*  -  2g 


(42) 


The  hydraulic  mean  depth  at  any  x-station  can 


.  _A  .  25SP  ,  _e_  r 

P  (8+1) 2g  6+1  ** 

w 


(43) 


where  £  is  the  corresponding  maximum  fluid  depth  at  that  station. 

Assuming  the  liquid  surface  is  open  to  the  atmosphere,  we  can  take 
p^  =  p^  ■  0  in  equation  (40).  Thus,  utilizing  equation  (43),  Bernoulli's 

equation  for  any  point,  x,  may  be  written: 


r_  +  i+i  „  i+i 

2«c  e  B 


g 

Dividing  equation  (44)  by  £*,  multiplying  by  — r 

A2 

number  as  Fr  *%/ — TZ  >  one  obtains: 
gcC* 


C*  (44) 

,  and  defining  the  Froude 


.  fo  .  g+l_jo.  ,  fo 

0+12  5*  ,  5 

6+1  ^ 


(45) 


Now  apply  the  energy  equation  to  a  compressible  fluid:  Heat  Added  +  Work 
Done  on  The  Fluid  *  Increase  of  Fluid  Energy.  Assuming  the  fluid  does  no 

P1  p2 

work,  the  work  done  on  the  fluid  (per  unit  mass)  is:  -  .  Neglecting 

P1  °2 

field  forces,  the  corresponding  increase  of  fluid  energy  per  unit  mass  is: 
12  12 

(e2  +  2  V2^  ~  ^el  +  ~2  Vi^ »  where  e^  ^  represents  internal  energy.  Finally, 

considering  an  adiabatic  process  (heat  added  -  0) ,  the  energy  equation  may 
be  written: 


r  i  Ft  1  0  2 

+  e.)  -  (-=■  +  e,)  -  y(v:  -  Vp 

1  P£  2  i  l  1 


(46) 
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and,  defining  enthalpy  H 


+  e. 


V1  V2  V2 

H! +  r  -  h2  +  r =  Ho  - H  +  r 


(47) 


Now,  for  a  perfect  gas,  (thermally  and  calorically) ,  H  *  c^T.  Furthermore, 
the  local  speed  of  sound,  c,  can  be  expressed  in  terms  of  temperature  T  by 


the  relation  c  ■  yFT,  where  R  *  c  -  c  is  the  gas  constant  for  air  and 

c  2  HR  H(c0-O  P  V 

y  *  -*•.  Thus,  c  ■  —  ■  - 2 - —  ■  H(y-l) ,  and  equation  (47)  becomes: 

c  c  c 

V  V  V 

2 

2  „2  c 

k  +  r  =  (A8) 

2  V 

Dividing  equation  (48)  by  c  and  defining  the  Mach  Number  as  M  *  — ,  one 

obtains: 


1  M2  _1_  YRTo 

Y-l  2  =  y-1  YRT  * 


or, 


(Y-l)  +  1 


(49) 


Note  that  equation  (49)  would  be  identical  to  equation  (45)  if  the  Mach 

T 

number  were  equal  to  the  Froude  number,  if  the  temperature  ratio  corres¬ 
ponded  to  the  depth  ratio  -jp-,  and  if  the  following  relationship  held: 


8+1 


(Y-l)  • 


(50) 


It  thus  follows  that  depth  changes  due  to  wave  motion  in  an  incompressible 
open  channel  flow  are  directly  proportional  to  temperature  changes  due  to 
wave  motion  in  a  compressible  gas-dynamic  flow;  and  that  in  the  hydraulic 
analogy  there  is  a  one-to-one  correspondence  between  the  incompressible 
Froude  number,  which  measures  the  importance  of  gravity  in  the  wave  motion, 
and  the  compressible  Mach  number,  which  measures  the  importance  of  fluid  com¬ 
pressibility  in  the  wave  motion.  To  complete  the  analogy,  we  introduce  the 
perfect  gas  relationship,  p  •  pRT,  to  obtain,  using  equation  (39), 
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sc  that  once  again  we  note  the  correspondence  between  pressure  waves  in  a  com¬ 
pressible  gas  and  gravity  waves  on  the  free  surface  of  a  liquid.  It  might  be 
well  at  this  point  to  review  the  assumptions  which  are  inherent  in  the  hy¬ 
draulic  analogy,  and  to  summarize  the  results: 

Assumptions :  a)  Bernoulli's  equation  (and  to  some  extent,  the  continuity 
equation),  as  used  here  for  an  incompressible  fluid,  pre-supposes  the  exis¬ 
tence  of  an  isotropic,  homogeneous,  inviscid,  Newtonian  continuum,  in  a  state 
of  thermodynamic  equilibrium,  undergoing  steady,  essentially,  one-dimensional 
flow  in  a  rigid  channel  of  uniform,  constant  cross-section,  where  gravity  is 
the  only  body  force,  and  where  there  is  no  heat  transfer,  no  significant 
change  in  internal  energy,  no  significant  energy  losses,  and  no  work  done  on 
the  fluid  (by  a  pump,  for  example)  or  by  the  fluid  (on  a  turbine  or  other 
device) .  b)  The  energy  equation,  as  used  here  for  a  thermally  and  calorically 
perfect  gas,  also  pre-supposes  the  existence  of  an  isotropic,  homogeneous,  in¬ 
viscid,  Newtonian  continuum,  in  a  state  of  thermodynamic  equilibrium,  under¬ 
going  steady,  essentially  two-dimensional  flow  through  a  region  of  uniform 
cross-section,  where  gravity  is  the  only  body  force,  and  where  there  is  no 
heat  transfer,  no  significant  energy  losses,  and  no  work  done  by  the  fluid, 
c)  The  rigid  channel  is  assumed  to  be  wide,  and  to  have  a  uniform,  symmetric 
constant  cross-section  which  can  be  defined  by  a  power-law  relationship  and 
which  is  open  to  the  atmosphere  at  its  top.  The  fluid  in  the  channel  is  fur¬ 
ther  assumed  to  maintain  a  shallow,  relatively  constant  (in  the  x-direction) 
depth  and  the  effects  of  surface  tension  are  neglected. 

Over  the  years,  many  investigators  have  addressed  themselves  to  the 
consequences  and  implications  of  this  rather  impressive  list  of  assumptions 
as  they  relate  to  the  validity  of  the  hydraulic  analogy  (see,  for  example, 
Preiswerk,  1940,  Ippen  and  Harleman,  1950,  and  Harleman  and  Crossley,  1952). 
Suffice  it  to  say  that  the  analogy  is  far  from  universally  valid,  especially 
in  situations  where  the  continuum  assumption  breaks  down  (rarefied  gas 
dynamics  in  the  upper  atmosphere)  or  where  there  are  substantial  losses  of 
energy  (such  as  when  shock  waves  form  in  supersonic  flows).  However,  even 
in  these  extreme  situations,  "the  water  table  has  been  strikingly  successful 
in  reproducing  most  qualitative  aspects  of  high  speed  flows.  It  has  the 
great  advantages  of  simplicity,  low  cost,  rapid  operation,  and  gives  the 
operator  a  chance  to  'play'  with  various  flow  patterns.  It  is  especially 
suitable  for  observing  unsteady  and  transient  flow  patterns."  (Shapiro, 

1949) .  One  might  add  that  the  water  channel  lends  itself  much  more  conven¬ 
iently  and  cheaply  to  effective  flow  visualization  techniques  than  does  any 
other  form  of  experimental  procedure.  Moreover,  in  the  case  of  shock  phenom¬ 
ena,  the  analogy  can  be  corrected  by  replacing  the  gas-dynamic  equations 
with  appropriate  shock  relationships  and  the  steady  flow  hydraulic  equations 
with  hydraulic  jump  equations. 

In  summary,  "Despite  the  many  questionable  features  of  the  analogue, 
the  water  channel  is  a  valuable  auxiliary  tool  which  has  the  ability  to  shed 
light  on  problems  of  current  importance  in  two-dimensional  flow.  It  justi¬ 
fies  itself  through  economy  of  time  and  effort  and  is  useful  for  making 
qualitative  studies  and  quantitative  experiments  of  a  comparative  nature" 
(Shapiro,  1949) . 
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With  this  in  mind,  table  1  has  been  prepared  to  highlight  the  essential 
features  of  the  analogy  which  allows  aerodynamic  situations  to  be  studied 
in  open  channel,  shallow  water  models. 

We  call  particular  attention  to  the  analogue  relationship,  equation  (50), 
which  establishes  the  cross-section  geometry  required  of  the  open  channel  as 
a  function  of  the  specific  heat  ratio  of  the  gas  to  be  studied.  Equation  (50) 
is  plotted  in  Figure  9,  where  it  is  observed  that  the  curve  is  asymptotic  to 
Y  *  2  for  8  -+■  <*>.  Rearranging  equation  (32)  reveals  that  8  -*■  •  corresponds  to 
the  condition: 


[ 


x(y)il/e 

5m 


£_ 

«m 


1  for  all  x  >  0,  or, 
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constant  for  all  x  bigger  than  zero.  If  (^— )  <  1,  then 
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(^— )®  -*■  0  as  8  00 ,  so  this  condition  only  holds  along  the  special  line  x(y) 

®m 

*  0.  Thus,  we  have,  for  x(y)  ■  0,  (^“)  <,  1  and  for  x(y)  *  arbitrary,  *  1, 

8m  8m 

which,  together,  define  a  channel  of  rectangular  cross-section,  having  half- 
width  and  arbitrary  depth  Note  also  that,  for  such  a  channel,  C  “  5*, 


as  observed  by  letting  8  -*■  00  in  equation  (43).  That  is,  the  hydraulic  mean 
depth  for  a  shallow  channel  of  rectangular  cross-section  is  simply  the  same 
as  the  actual  fluid  depth  at  any  x-station. 


In  most  cases,  the  use  of  rectangular  channels  is  convenient,  but,  of 
course,  assumes  the  existence  of  a  hypothetical  perfect  gas  having  a  specific- 
heat  ratio,  y,  equal  to  2.  For  most  real  gases,  y  is  more  nearly  in  the 
range  1.2  to  1.5  (Streeter,  1961).  Its  value  in  air,  for  example, 
is  1.4,  which  corresponds  to  a  channel  having  a  cross-sectional  contour  given 
by: 


„  .  Xkl  .  (X-,2/3  .  (52, 

m  8m 

Such  a  channel  is  illustrated  in  Figure  10,  as  are  the  corresponding 

cross-sections  for  y  »  1.5  (triangular  contour,  8  ■*  1) ,  y  *  1*2  (8  *  ^)  and 

Y  ■  2  (rectangular  cross-section) .  Note  that  the  triangular  contour,  which 
is  reasonably  convenient  to  design  and  build,  provides  a  good  approximation 
(y  *  1.5)  to  the  analogy  for  aerodynamics  (y  ■  1.4).  Fortunately,  however, 
it  turns  out  that  the  rectangular  channel  (y  “  2.0),  which  is  even  more 
convenient  from  the  point  of  view  of  scaling  the  body  models  to  be  studied 
(Whicker,  1951) ,  gives  results  that  are  quite  acceptable  for  qualitative 
flow  assessment.  In  fact,  the  errors  thus  introduced  are  not  too  serious  in 
both  the  subsonic  and  in  the  supersonic  regimes  (Streeter,  1961)  and  they 
are  less  important  than,  for  example,  the  boundary  layer  which  develops 
on  the  bottom  of  the  fluid  stream,  the  viscosity,  the  thermal  conductivity 
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Plotted  as  a  Function  of  Perfect  Gas  Specific  Heat 
Ratio,  v. 
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Figure  10:  Open  Channel  Cross-Sectional  Geometries  for  the  Perfect 
Gas  Specific  Heat  Ratios  1.0,  1.2,  1.4,  1.5  and  2.0. 


and  the  surface  tension.  For  this  reason,  it  was  decided  to  pursue  the 
experimental  aspects  of  this  investigation  utilizing  flow  visualization 
and  hydraulic  analogue  measurements  in  rectangular  water  channels  as 
described  in  the  section  that  follows. 


SECTION  VI 

EXPERIMENTAL  STUDIES  OF  WAKE 
DEVELOPMENT  DOWNSTREAM  OF  A 
DOUBLE-CYLINDER  CONFIGURATION 


Experimental  Arrangement 


For  this  phase  of  the  investigation,  a  TECHNOVATE  Hydrodynamics  System, 
Model  9093  was  employed  (Scott-Engineering  Sciences,  a  division  of  A-T-0, 

Inc.,  910  Southwest  12th  Avenue,  Pompano  Beach,  Florida  33060).  This  system, 
shown  in  Figure  11,  consists  of  a  shallow  rectangular  channel  in  which  water 
is  recirculated  by  a  centrifugal  pump,  a  precision  level  gage  for  wave  height 
measurements,  and  a  hydrogen  bubble  generator  to  allow  for  flow  visualization. 
The  channel  working  section  measures  18  Inches  (45  3/4  cm.)  wide  by  2  1/2 
inches  (6  1/3  cm.)  deep  by  48  inches  (122  cm.)  long.  Thus,  the  maximum  depth- 
to-width  ratio  is  0.139  which  reasonably  satisfies  the  shallow,  wide  channel 
requirements  of  the  hydraulic  analogy.  Moreover,  the  inlet  tank  and  channel 
entrance  are  carefully  profiled  to  provide  satisfactory  flow  conditions,  as 
are  the  working  section  length  and  discharge  configurations.  All  are 
supplied  in  the  form  of  an  accurate  reinforced  glass  fiber  moulding  which 
ensues  a  uniform  smooth  cross-section  throughout. 

The  flow  system  is  mounted  by  levelling  screws  on  a  robust  fabricated 
steel  frame  fitted  with  casters  for  maneuverability.  Beneath  the  frame,  a 
centrifugal  pump,  controlled  by  a  diaphragm  type  valve  drives  the  channel 
fluid  in  a  continuously  recirculating  pattern.  The  maximum  pump  output 
delivers  18  U.S.  gallons  per  minute  which  is  equivalent  to  a  volumetric  flow 
rate  of  .04  ft^/sec  (7.48  U.S.  gallons  ■  1  cubic  foot)  or  1135.7  cm^/sec. 

The  minimum  pump  flow  can  be  throttled  all  the  way  down  to  zero  by  controlling 
the  diaphragm  value.  Thus  volumetric  flow,  $,  through  the  working  section 
is  given  by: 

0  <_  <_  1135.7  cm^/sec 

A 

which  corresponds  to  an  average  fluid  velocity  of  7-,  where  A.  is  the  fluid- 

1  1 

filled  cross-sectional  area  of  the  working  section.  Note  that  A^  -  2gj£*, 
where  2g^  *  the  channel  width  •  45.75  cm,  and  £*  ■  the  mean  depth  of  the 
flowing  liquid.  Now  recall  that  the  Froude  Number  is  defined  by: 


Fr 


ir 
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■J. 
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4gJc*Vgc 


Therefore,  Che  relationship  becween  Froude  Number,  volumetric  flow  and  fluid 
depth  for  this  hydrodynamics  system  is  (g£  *  981.5  cm/sec^): 


Fr 


5*3(2.05  x  106) 


(53) 


3  (- 3/21 

For  4*0,  Fr  ■  0  and  for  4  •  4  m  1135.7  cm  /sec,  Fr  “  .7935*  •  In 

max 


general,  for  a  desired  Froude  number,  the  fluid  depth  may  be  varied  depending 
upon  the  volumetric  flow  rate  at  which  one  chooses  to  operate  the  system.  How¬ 
ever,  this  is  not  entirely  arbitrary,  since  the  operating  fluid  depth  must 
further  take  into  account  the  growth  of  a  boundary  layer  along  the  bottom  of 
the  channel  —  a  viscous  effect  which  is  neglected  in  the  hydraulic  analogy. 


Hinze  (1959  pg.  487)  has  suggested  that  the  following  equation  may  be 
used  to  obtain  an  estimate  of  the  increase  in  boundary  layer  thickness,  6, 
with  downstream  distance,  x,  from  the  origin  of  a  flat  surface: 
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0.045 
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0.79 


(54) 


where  6  «  0.126  defines  the  momentum  thickness  and  v  is  the  kinematic 

m 

viscosity  of  the  fluid.  Rearranging  equation  (54)  and  replacing  5^  by  5 
yields: 
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0.375 


(55) 


Asstiming  the  fluid  is  water  at  room  temperature  (68°F  *  20°C) ,  its  kinematic 
viscosity,  v,  is  equal  to  1.004  centistokes  (Bolz  and  Tuve,  1970,  pg.  67) 
which  is  .01004  cm^/sec.  Furthermore,  with 
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2gL5* 
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45.755*’ 


equation  (55)  becomes: 
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0.318 


30.21 


(56) 


But,  from  the  earlier  definition  of  the  Froude  number,  i.e.,  equation  (53): 


1432.33  (Fr)  5*3/2  -  4 


(57) 
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From  here,  we  see  that,  for  a  given  depth,  £*,  the  Froude  nuober  increases 
directly  with  flow  rate  4.  Thus,  at  any  operating  depth,  the  maximum  Froude 
number  which  can  be  achieved  is  limited  by  the  maximum  flow  (1135.7  cm -V sec) 
that  the  recirculating  pump  can  deliver.  Substituting  this  value  of  ♦  into 
equation  (57)  yields  equation  (58)  which  is  plotted  in  Figure  12: 

(Fr)  -  . 793(£*)~3/2  (58) 

max# 


Equation  (58) ,  or  Figure  12  both  Indicate  that  to  get  up  to  a  Froude  number 
of  at  least  unity  in  this  flow  system,  the  fluid  depth  must  be  reduced  to 
about  8  millimeters  at  maximum  pump  output.  It  we  incorporate  this  infor¬ 
mation  into  equation  (56),  i.e.,  let  $  -  1135.7  cm^/sec  and  £*  *  .857  cm, 
we  find  that  the  growing  boundary  layer  thickness,  6,  has  already  reached 
fluid  depth  £*  when  x  -  23.7  cm,  or,  about  202  of  the  way  into  the  channel. 
Beyond  this  point,  the  entire  flow  is  in  the  boundary  layer;  and  the  situation 
worsens  at  lower  Froude  numbers  because  as  $  is  reduced  for  £*  held  constant, 
the  boundary  layer  thickens  sooner  with  increasing  x.  The  dimensions  and 
flow  characteristics  of  this  particular  water  table  thus  seem  to  Introduce  a 
problem  in  terms  of  one  of  the  assumptions  Inherent  in  the  hydraulic  analogy. 
To  examine  further  the  significance  of  this  problem,  let  $  ■  1135.7  cm3/sec 
and  combine  equations  (58)  and  (56)  to  get: 
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6  «  0.318{^-}‘79  (121.92)  *79  5*  [1x35.7]  *21 
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Furthermore,  writing:  (£*)”'79  *  [(£*)”^’^] *^27 

Fr  ,„7 

*  (  793^ "  from  equation  (58) , 

therefore:  -  3.647{r-}0’79  (Fr  ) ’527  (59) 
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where  x^  *  121.92  represents  the  total  length  of  the  channel.  Again,  note 


that  when  (Fr)  *  1,  6/C*  ■  1  for  —  * 
max  x. 

6  X 

the  other  hand,  if  we  set  -rr  ■  1  for  — 

*L 


0.194  as  already  mentioned.  On 

1  1,  then  we  find  that  (Fr)  * 

max 


.086,  which  is  much  tod  low  for  our  purposes.  Any  way  that  one  views  this 
situation,  it  clearly  represents  a  source  of  error  in  the  experimental 
arrangement.  The  extent  of  this  error  in  terms  of  the  hydraulic  analogy 
will  require  further  investigation,  but  its  magnitude  in  the  present  flow 
system  is  highlighted  by  equation  (59),  which  is  also  plotted  in  Figure  13 
with  (Fr)  acting  as  a  parameter.  From  these  results,  one  would  have  to 
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conclude  that  all  measurements  made  in  the  channel  working  section  are 
being  taken  within  the  boundary  layer  of  the  flow.  One  possible  way  around 


Maximum  Froude  Number  vs.  Water  Depth  for  the 
Model  9093  Water  Table. 
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BOUNDARY  LAYER  GROWTH 


Boundary  Layer  Growth  Along  the  Bottom  of  a  Wide,  Shallow,  Open  Rectangular 
Channel,  as  a  function  of  Froude  Number;  (for  Model  9093  Test  Section). 


this  dilemma  is  to  reduce  the  width  of  the  channel  using  mountable  plexi¬ 
glass  plates  supplied  with  model  9093.  This  would  allow  a  deeper  flow  for 
the  same  cross-sectional  area,  but  might  adversely  affect  the  depth-to-width 
ratio,  which  must  be  kept  small.  It  is  thus  an  option  that  must  be  con¬ 
sidered  carefully.  Two  other  options  are  described  below. 

The  growth  of  a  boundary  layer  arises  from  frictional  interactions 
between  a  moving  fluid  and  a  solid  surface.  One  consequence  of  such  inter¬ 
actions  is  a  generalized  tendency  for  the  flow  in  the  test  section  to  be 
decelerated.  This  effect,  however,  may  be  cancelled  approximately  by  intro¬ 
ducing  a  gravity  force  obtained  through  tilting  the  channel  by  a  small 
amount.  The  proper  amount  of  tilt  can  be  found  most  practically  by  trial 
and  error.  However,  it,  too,  must  be  kept  small  because,  in  arriving  at 
the  analogue  relationships,  it  was  assumed  that  the  water  flow  is  two-dim¬ 
ensional.  This  assumption  rules  out  significant  vertical  components  of 
velocity  and  changes  in  water  depth  (the  latter  suggesting  the  use  of  large 
models  compared  with  the  water  depth) . 

As  a  final  point  on  boundary  layer  effects,  it  might  be  mentioned  that 
viscous  phenomena  such  as  boundary  layer  growth  may  be  accounted  for,  at 
least  approximately,  by  matching  as  closely  as  possible  the  Reynolds  Number 
of  the  Water  Table  Model  Tests  with  that  of  the  Aerodynamic  Prototype  experi¬ 
ments.  Though  this  may  often  be  difficult  or  impractical,  one  must  still 
remember  the  point  stressed  earlier  -  that  is,  that  a  water  channel  can  never 
give  precise  information  on  the  behavior  of  a  compressible  fluid,  but  is  best 
confined  to  a  qualitative  demonstration  of  the  flow  in  or  around  a  given 
model.  Its  major  asset  is  that  it  can  be  used  cheaply  and  conveniently  to 
assess  the  relative  merits  of  several  versions  of  a  design,  prior  to 
carrying  out  more  exact  investigations  in  a  sub-or-supersonic  wind  tunnel. 


Flow  Visualization  and  Surface  Wave  Measurements 


In  the  model  9093  Hydrodynamics  System,  the  flow  can  be  visualized 
either  by  the  Hydrogen  Bubble  Technique  (Schraub,  Kline,  et.  al. ,  1965  and 
Davies  and  Fox,  1967)  or  by  photographing  surface  wave  patterns  (Shapiro, 
1949) .  The  former  utilizes  as  flow  tracers  small  hydrogen  bubbles  evolved 
in  water  by  electrolysis  at  a  fine  wire  cathode  which  is  positioned  normal 
to  the  fluid  flow  (Figure  14).  By  careful  illumination  and  photography,  the 
paths  traced  out  by  these  bubbles  as  they  are  carried  along  by  the  flow  can 
be  observed.  Owing  to  their  small  size  (.0005  to  .0015  -  inch  diameter) 
they  follow  the  flow  accurately  and  their  rate  of  rise  due  to  the  natural 
buoyancy  of  the  light  gas  is  small  compared  with  the  average  rate  of  flow 
observed.  Hence,  the  flow  visualization  obtained  is  sufficient  for  many 
purposes  including  those  of  this  investigation. 

The  cathode  consists  of  a  fine  stainless  steel  wire  (.002-inch  diameter) 
supported  in  tension  by  a  two-pronged  fork  holder.  The  holder  is  grasped 
by  a  crocodile  clip  attached  to  the  end  of  a  rod  which  completes  the 
electrical  connection  to  the  wire.  The  electrical  lead  hooked  to  the  rod 
originates  from  the  negative  terminal  of  the  pulse  generator.  This  entire 
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Figure  14:  Hydrogen- Bubble  Generating 
Unit  on  Hydrodynamics 
Model  9093  Water  Table 


cathode  arrangement  can  be  positioned  where  required  in  the  channel  by  a 
combination  of  clamps  and  rod  holders  attached  to  the  tripod  shown  in 
Figure  14.  The  anode  is  located  at  the  channel  discharge  end. 

Since  ordinary  tap  water  can  have  a  wide  range  of  conductivities,  the 
addition  of  a  neutral  salt  is  recommended  in  order  to  ensure  the  generation 
of  a  large  number  of  small,  equal-sized  bubbles  which  are  uniformly  pro¬ 
duced  at  various  points  along  the  wire.  For  this  purpose,  three  liters  of 
a  saturated  solution  of  sodium  sulfate  is  added  to  the  fluid  in  the  channel, 
yielding  an  overall  salt  concentration  of  25  grams  per  liter  of  channel 
fluid. 

The  hydrogen  bubble  generator  produces  "on"  and  "off"  electrical  pulses 
with  provision  for  varying  the  off  period,  controlling  the  on  period  (pulse 
duration),  and  adjusting  the  D.C.  Voltage.  The  adjustable  time  base  permits 
quantitative  velocity  data  analysis  to  be  obtained  by  photographic  means; 
while  the  variable  D.C.  Voltage  is  used  to  control  bubble  size,  which  is 
roughly  propotional  to  current  density.  Normally,  current  densities  of  1  to 
2  amp/in^  (6  to  12  ma  per  inch  of  wire)  are  sufficient  to  produce  .001  to 
•0015-inch  diameter  bubbles  from  a  .002-inch  diameter  cathode  wire  -  and 
this  is  relatively  independent  of  the  velocity  of  flow  past  the  wire. 

To  illuminate  properly  the  hydrogen  bubbles  for  photographic  purposes, 
it  is  essential  to  provide  a  considerable  amount  of  light  in  a  fairly 
narrow  flat  beam,  roughly  parallel  to  the  water  surface  (and  below  it)  in  a 
plane  subtending  between  90°  and  110°  to  the  line  of  sight  of  the  observer. 
This  is  achieved  in  model  9093  by  providing  a  self-contained  light  source 
(a  cooled  55  watt  tungsten  iodine  bulb  backed  by  a  concave  mirror)  and 
light  guide  (made  of  polished  plexiglass)  which  enables  the  user  to  illum¬ 
inate  almost  any  part  of  the  flow  channel  using  the  principle  of  total 
internal  reflection.  An  example  of  what  can  be  achieved  is  shown  in  Figure 
15.  Here,  two  stainless  steel,  1"  diameter  cylinders  have  been  immersed  in 
the  water  channel  to  simulate  the  condition  illustrated  in  Figure  1.  The 
cylinders  are  rigidly  fastened  together  as  the  flow  is  accelerated  to  show 
the  double-vortex  formation  in  the  wake.  Photographic  records  were 
obtained  using  a  400  ASA  film  with  a  wide  aperture  lens  and  an  exposure  of 
1/30  th  at  f  2.8. 

The  flow  visualization  technique  involving  photography  of  surface  wave 
patterns  is  based  upon  the  fact  that  when  diffuse  light  is  reflected  off 
of  the  flow  surface,  the  troughs  of  the  waves  appear  much  brighter  than  the 
crests.  Since  the  regions  with  higher  waves  correspond  in  the  analogy  to 
regions  of  high  pressure,  and  vice  versa,  the  light  gradations  on  photo¬ 
graphic  records  of  the  wave  patterns  can  be  calibrated  to  yield  information 
about  the  pressure  distribution  in  the  flow  field.  Moreover,  this  technique 
also  provides  useful  qualitative  descriptions  of  the  formation  and  behavior 
of  streamline  patterns  and  shock  waves  around  different  geometric  config¬ 
urations  which  intercept  the  flow. 

Photography  of  wave  patterns  in  the  model  9093  can  be  accomplished 
conveniently  by  using  high  intensity  spot  lights;  in  this  case,  Berkey 
Colortran  model  112-031  cine-queen  1000  watt  lamps.  A  demonstration  of  the 
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Figure  15:  Hydrogen  Bubble  Flow  Visualization  Downstream  of  Two, 

Tangent  Stainless  Steel,  One- Inch  Diameter  Cylinders 

Immersed  in  the  Water  Channel  with  the  Flow  Travelling 

From  Left  to  Right. 

i 

\  4, 


types  of  results  that  can  be  obtained  with  this  lighting  system  is  shown  in 
Figure  16.  Here,  the  spot  lights  have  been  reflected  off  the  ceiling  of 
the  laboratory  to  produce  a  diffuse  light  and  thus  eliminate  glare  from  the 
water  surface.  Moreover,  the  lights  have  been  positioned  so  that  the  dif¬ 
fuse  beam  is  reflected  from  the  water  surface  directly  into  the  camera  lens 
when  the  water  is  absolutely  still.  This  gives  the  observer  a  reference  to 
the  angle  or  slope  of  the  water  surface  at  any  particular  point.  Recall 
that  the  subsequent  troughs  of  the  waves  appear  to  be  much  brighter  than  the 
crests. 

The  flow  photographed  in  Figure  16  is  interacting  with  a  simple  plexi¬ 
glass  model  of  the  arm  and  torso  configuration.  This  model,  shown  in  Figure 
17,  was  scaled  and  constructed  from  anthropometric  data  such  that  the  ratio 
of  its  frontal  cross  sectional  area  to  that  of  the  flow  channel  was  17:1. 

The  arms  were  attached  with  wire  loops  passed  through  the  eye-hooks,  so  that 
they  could  "swing"  freely  about  the  shoulders.  The  only  restraint  imposed 
was  in  the  form  of  a  plexiglass  back  rest  that  was  built  into  the  apparatus 
which  held  the  model  in  place  in  the  water  table. 

Figure  16  has  been  extracted  for  illustrative  purposes  from  a  20-minute 
video-tape  that  was  prepared  to  show  the  flow/model  interactions  for  a  range 
of  Mach  numbers.  The  Mach  number  in  this  particular  photograph  is  0.85,  and 
several  features  should  be  noted.  First,  the  large  gap  between  wave  patterns 
indicates  a  significant  transition  from  a  high  pressure  region  in  the  front 
of  the  model  to  a  low  pressure  region  behind  it.  Second,  one  can  see  a 
developing  bow  shock  wave,  with  a  corresponding  increase  in  pressure  over 
the  forward- facing  parts  of  the  model.  Additional  developing  shocks  (across 
which  there  is  a  fall  in  pressure)  are  in  evidence  along  the  wider  parts  of 
the  model  and  at  the  trailing  edges,  where  the  pressure  rises  again.  In 
particular,  observe  how  the  "arms"  in  this  model  have  been  thrown  clear  of 
the  torso,  and,  although  qualitative,  how  clearly  one  can  visualize  the 
developing  flow  patterns  as  an  obstruction  is  placed  in  the  oncoming  flow. 
This  is  one  of  the  key  features  of  this  technique  as  a  means  for  under¬ 
standing  fluid  behavior  in  a  complex  situation.  Quantitative  information 
concerning  actual  water  level  (hence,  corresponding  fluid  pressure)  cam  be 
further  obtained  by  means  of  a  built-in  precision  level  gage  which  is  an 
integral  feature  of  model  9093.  The  gage  has  a  calibrated  depth  probe  which 
causes  an  indicator  light  on  the  control  panel  to  come  on  as  soon  as  the 
probe  tip  contacts  the  water  surface.  More  details  of  the  results  observed 
on  the  20-minute  video  tape  are  discussed  below. 


Results 


In  the  experiments  to  be  discussed  the  model  shown  in  Figure  17  was 
held  in  place  by  a  series  of  three  laboratory  clamps.  These,  in  turn,  were 
connected  to  the  supportive  structure  which  runs  along  the  runners  on  the 
channel  (see  Figure  14).  A  five  pound  counter-weight  was  used  to  hold  the 
model  above  the  water  such  that  it  could  be  launched  into  a  fully  developed 
flow.  The  flow  was  controlled  by  the  diaphragm  valve  described  previously 
and  the  actual  water  speed  was  measured  by  timing  a  free-floating  particle 
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as  it  traversed  through  a  known  channel  distance.  Two  separate  series  of 
experiments  were  performed:  in  the  first,  the  model  was  positioned  in  the 
water  and  the  speed  of  the  flow  was  increased  gradually  from  zero  to  one 
which  corresponds  to  an  aerodynamic  flow  of  Mach  1;  in  the  second,  the 
model  was  actually  launched  into  pre-set  and  fully  developed  flows  corres¬ 
ponding  to  the  same  set  of  Mach  numbers,  i.e.,  0  to  1.  Results  of  the  first 
set  of  experiments  are  described  in  the  next  paragraph. 

As  the  flow  was  gradually  increased  from  zero,  wave  patterns  formed 
first  in  front  of  the  model,  where  stagnation  points  developed  as  fluid 
momemtum  was  shifted  90°.  The  force  build-up  at  such  stagnation  points  had 
a  tendency  to  pull  the  arms  away  from  the  torso,  and  this  was  clearly  illus¬ 
trated  in  the  video  tape.  As  the  liquid  flow  rate  continued  to  increase, 
the  water  level  (£)  continued  to  build  up  against  the  front  of  the  model, 
indicating  correspondingly  higher  pressures  (see  equation  [51])  in  the 
analogous  aerodynamic  system.  Behind  the  model,  in  the  wake  region,  the 
water  level,  and  hence  fluid  pressures  were  much  lower  (see,  for  example. 
Figure  15  and  16).  Eventually,  the  model  became  totally  submerged  and  some 
vibration  of  the  arms  was  noted.  Vibration,  or  at  least  some  slight  oscil¬ 
lation  in  the  arms  was  first  observed  at  Mach  0.4.  Recall  that  these  are 
caused  in  part  by  the  small  opening  created  when  the  upper  arms  pull  away 
from  the  torso,  thus  generating  a  Venturi- like  effect  which,  coupled  with 
gravity  and  the  vortex  street  generated  by  this  blunt-body  configuration, 
drives  the  arms  back  towards  the  torso  . . .  where  stagnation  forces  cause  the 
process  to  repeat  itself.  At  Mach  0.6,  the  arms  no  longer  touched  the  torso 
but  still  continued  to  vibrate  slightly.  The  possibility  of  resonance 
effects  coming  into  play  due  to  phenomena  such  as  these  was  suggested  as 
possibly  being  influential  in  causing  the  ejection  seat  occupant  to  lose 
control  of  his  or  her  arms,  and  flail  injury  may  very  will  be  the  immediate 
result.  At  Mach  0.8,  the  limb  displacement  seemed  to  decline  somewhat  and 
the  vibrations  occurred  at  a  lower  frequency.  This  may  have  been  due  to  the 
increasing  frictional  interaction  between  the  arms  and  the  back  seat  rest, 
and  to  less  influence  from  gravitational  effects  because  of  the  time  scales 
involved.  When  the  flow  rate  reached  Mach  1,  the  model  became  completely 
submerged  and  virtually  invisible  to  the  above  water  observer.  The  water 
height  (?  =  3  to  4  times  5Q)  clearly  indicated  the  presence  of  large  forces 
tending  to  drive  the  arms  away  from  the  torso,  and  indeed,  the  arms  flailed 
wildly  even  in  this  simple  model  of  the  actual  ejection  process.  Moreover, 
the  large  increase  in  water  level  in  front  of  the  model  was  accompanied  by  a 
similar  decrease  in  the  water  level  behind  the  model,  so  that  the  total 
pressure  drop  across  this  geometric  configuration  was  rather  substantial  for 
Mach  Numbers  in  excess  of  0.70.  There  was  no  arm-torso  contact  beyond  Mach 
0.5  and  a  continuous  stream  of  fluid  could  be  seen  flowing  in  between  the 
arms  and  the  torso  for  this  and  higher  Mach  numbers. 

Turning  our  attention  now  to  results  obtained  utilizing  the  "launching" 
system  described  previously,  it  was  observed  that  a  simulated  ejection 
taking  place  at  Mach  0.4  has  associated  with  it  the  development  of  some 
rather  significant  high  pressure  regions  fore,  and  low  pressure  regions  aft 
of  the  pilot-ejection-seat  configuration.  The  pressure  drop  across  the 
configuration  was  estimated  to  be  some  four  times  that  of  the  free  stream 
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dynamic  pressure  (i.e.,  C  =  -50)»  which  compares  favorably  with  values 
predicted  by  the  theory  developed  earlier.  Furthermore,  during  this  series 
of  experiments,  it  was  again  observed  that  events  start  to  get  critical  when 
ejection  Mach  numbers  exceed  0.60.  At  this  speed,  the  arms  were  thrown  out 
violently  when  the  model  was  launched  suddenly  into  the  flow  and  some 
slight  oscillation  of  the  limbs  was  noted.  By  the  time  the  Mach  number 
reached  0.8,  a  clearly  defined  oscillating  wave  pattern  was  set  up  behind 
the  model  due  to  the  vortex  street  generated  by  flow  separation  effects. 

Also,  the  early  stages  of  the  formation  of  a  bow  shock  wave  could  be  seen  as 
the  fluid  piled  up  in  front  of  the  model  and  successive  crests  of  the  bow 
waves  became  compressed,  i.e.,  appeared  closer  together  than  before.  Moving 
towards  the  rear  of  the  model,  one  could  see  large  gaps  between  wave  patterns, 
indicating  a  rather  significant  transition  from  a  high  pressure  region  in 
front  to  a  very  low  pressure  region  behind  the  model.  Quantitative  informa¬ 
tion  concerning  actual  water  level  was  obtained  on  the  model  9093  by  means 
of  the  built-in  precision  level  gage  and  also  by  means  of  a  specially 
designed  small  wave  height  measuring  device  (Bishop,  G.  H. ,  1980  and  Cahill, 

B. ,  1981) .  Both  methods  revealed  the  presence  of  maximum  pressure  gradients 
on  the  order  of  5  to  8  times  that  of  the  measured  free- stream  dynamic 
pressure,  again  agreeing  rather  well  with  theoretical  predictions  and 
verifying  that  the  pilot  ejected  under  similar  circumstances  would  most 
likely  suffer  serious  windblast  and/or  flail  injury.  At  Mach  1,  events  were 
still  further  magnified,  with  large  limb  displacements,  lower  oscillation 
frequencies  and  violent  shock/boundary- layer  type  interactions  that  were 
suggestive  of  serious  potential  consequences. 

In  summary,  the  experimental  segment  of  this  research  program  appeared 
to  corroborate  the  theoretical  formulation  that  predicted  significant  prob¬ 
ability  of  flail  for  ejections  occurring  above  Mach  0.70.  Based  on  these 
confirming  observations,  some  concluding  points  are  presented  in  the  next 
section. 


SECTION  VII 
CONCLUDING  REMARKS 


This  research  program  has  examined  the  serious  problem  of  windblast 
from  both  the  theoretical  and  experimental  points  of  view.  The  mathemat¬ 
ical  theory  developed,  although  greatly  simplified  for  purposes  of  tracta- 
bility,  nevertheless  revealed  the  essential  physics  which  leads  to  limb 
flailing.  That  is,  the  theory  clearly  identified  fluid  dynamic  stagnation 
and  flow  separation  forces  as  first-order  effects  which  contribute  to  the 
dislodgment  and  subsequent  flailing  of  the  limbs  of  an  ejection  seat  oc¬ 
cupant.  These  effects  are  strongly  dependent  upon  the  angle  at  which  the 
pilot  intercepts  the  air  stream,  and  they  become  significant  and  quite 
serious  for  ejections  taking  place  above  Mach  0.7.  To  confirm  these 
findings,  a  series  of  water  table  experiments  was  performed  in  which  flow 
stagnation  and  separation  effects  were,  indeed,  seen  to  be  primary  causes 
of  windblast  and  flail  phenomena.  Both  qualitative  and  quantitative 
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measurements  of  the  flow  patterns  revealed  that  the  critical  Mach  number  at 
which  the  aerodynamics  of  windblast  becomes  potentially  harmful  is  around 
0.6,  and  beyond  that  point  flail  is  almost  assured.  Parameters  of  signifi¬ 
cance  for  design  purposes  include  the  trajectory  of  the  ejection  process, 
the  geometry  of  the  ejection  configuration  and  the  speed  of  ejection.  All 
of  these  are  not  suited  for  supersonic  ejections  as  presently  designed  and 
significant  changes  appear  to  be  necessary  in  order  to  insure  the  safety  of 
the  aircrew  during  combat  missions. 
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